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ABSTRACT 

Peculiar velocities in the nearby Universe can be measured via the kinetic 
Sunyaev-Zel'dovich (kSZ) effect. Using a statistical method based on an opti- 
mised cross-correlation with nearby galaxies, we extract the kSZ signal generated 
by ionised coronae of galaxies from the Cosmic Microwave Background (CMB) 
temperature anisotropics observed by the Wilkinson Microwave Anisotropy 
Probe (WMAP). Marginalising over the thermal Sunyaev-Zel'dovich contribu- 
tion from clusters of galaxies, possible unresolved point source contamination, 
and Galactic foregrounds (dust, synchrotron and free-free emission), we report 
a kSZ bulk flow signal present at the ~ 90% confidence level in the seven-year 
WMAP data. When only galaxies within 50h^^ Mpc are included in the kSZ 
template we find a bulk flow in the CMB frame of \V\ — 533 ± 263 km s^^, in 
the direction / = 324 ± 27, b = — 7± 17, consistent with bulk flow measurements 
on a similar scale using classical distance indicators. We show how this compari- 
son constrains the (ionised) baryonic budget in the local universe. On very large 
(^ 500 h^^ Mpc) scales, we flnd a 95% upper limit of 470 km s^^, inconsistent 
with some analyses of bulk flow of clusters from the kSZ. We estimate that the 
significance of the bulk flow signal may increase to 3-5(7 using data from the 
Planck probe. 



1 INTRODUCTION 



Peculiar velocities are the only probe of the large-scale 
(10 - 1000 h-^ MpcQ distribution of mass in the nearby 
Universe. Surveys of the peculiar velocity field have re- 
cently returned to spotlight, in part due to the availability 
of larger surveys of galaxy distances, such as the SF1++ 
( Masters et al.pOOG i sample of TuUy-Fisher distance or 
the "First Amendment" supernova compilation (Turn- 
|bull et al.|20T2l ), as well as new analysis methods such as 
the minimum variance weighting method for bulk flows 



Background (CMB). Therefore, as the kSZ errors do not 
depend on distance, it is a promising probe of peculiar 
velocities on large scales. 

Some recent peculiar velocity studies using classical 
distance indicators have suggested that there may be a 
cosmic bulk flow on scales of ~ 50 — 100 Mpc in 



excess of that expected in ACDM (Watkins et al. 112009 



Lavaux et al. 2010 Feldman et al. 2010 Colin et al. 2011 1, 



whereas other studies have found lower values (Nusser 



Davis 2011 Turnbufl et al. 20121. Clearly, it would 



(Watkins et al. 2009 1 or the Monge-Ampere-Kantorovitch 
reconstruction ( [Lavaux et al.|[2bl0[ ). However, peculiar 
velocity surveys based on galaxy distances probe scales 
smaller than ~ 100/i~^ Mpc. The limitation is due to 
the linear increase of the peculiar velocity error with dis- 
tance, and the cost of collecting the large samples that 
are necessary to statistically reduce this error. In con- 
trast, the kinetic Sunyaev-Zel'dovich effect ( jSunyaev fc| 
Zeldovich|1980( hereafter kSZ) is sensitive to relative mo- 
tion of ionised gas with respect to the Cosmic Microwave 



^ The parameter h is today's Hubble constant in units of 
100 km s-l Mpc-^ 



be useful to have an independent probe of the peculiar 
velocity field on these scales. 

Still more controversially, on much larger (~ 
S00h~^ Mpc) scales, the kSZ effect in rich clusters has re- 
cently been used to measure the bulk flow by |Kashlinsky| 
hereafter KAKE). However subsequent stud- 



et al. 



(2008 



ies using similar data sets (Keisler 2009 Osborne et al. 
2011 Mody fc Hajian|2012 1 have found larger uncertain- 



ties than those claimed by KAKE. The main challenge is 
to disentangle the kSZ signal from the primary fluctua- 
tions of the CMB and the instrumental noise. The kSZ 
signal is at the level of a few fiK, while the CMB is in 
tens of /iK and the noise may be in hundreds of for 
WMAP ( |Jarosik et al.|2011| |. Previous kSZ work has fo- 
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cused on clusters of galaxies ( Haehnelt fc Tegmark|1996 
Kashlinsky et al.|2008||Osborne et al.|2011[|Mody fc Ha-| 



jian 2012 I because clusters hold a significant fraction (but 



not most) of the ionised plasma in the local universe. In- 
deed recently [Hand et al.| ( |201^ have detected the signal 
of pairwise infall of clusters using the kSZ effect. Unfortu- 
nately there are relatively few clusters of galaxies within 
100-200/i~^ Mpc in the nearby Universe, and so, as shown 
by the studies cited above, the corresponding errors on 
the bulk flow are large. 

In this paper, we use galaxies themselves to produce 
a kSZ template which we cross-correlate with the CMB 
fluctuations. Galaxies, including spirals, are expected to 



have a large plasma corona (Fukugita & Peebles 2006 
Rasmussen et al.|2009 \ , indeed this is a standard assump- 



tion of galaxy formation models for over 40 years (Rees| 
irOstriker,1977, Silk..l977, White fc Rees,1978J . While 
much of the hot gas is expected to lie outside halos, in the 



so-called warm-hot intergalactic medium (WHIM, Dave 
et al.||2001| iShuU et al.||2011l), nevertheless, this WHIM 



gas should be correlated with the haloes of larger galax- 
ies and thus may be detected statistically through cross- 
correlations (see Soltan|2006 for a claimed detection). 

The use of galaxies as kSZ tracers has also been ad- 



vocated by Ho et al.^(2009) and ^Shao et al. ( 2011 1. In par- 
ticular, using 2MASS galaxies as a probe of the KAKE 



bulk flow has been proposed by Zhang ( 2010 1. On the one 



hand, in a given volume there are far more galaxies than 
there are clusters of galaxies, e.g., within 200/i~^ Mpc, 
there are ~60 000 galaxies in the 2M-I— I- galaxy catalogue 



(Lavaux & Hudson 20111, whereas in the Reflex-eBCS- 
CIZA catalogue (RBC, [Kocevski fc Ebeling|2006[ ), there 
are only 273 clusters of galaxies in the same volume. On 
the other hand, the density profile of the plasma corona 
is far more uncertain for galaxies than for the clusters of 
galaxies. Moreover, the corona occupies a much smaller 
area on the sky, which requires an accurate modelling 
of beaming effects. Nevertheless there are expected to be 
more ionised electrons in or near galaxies than in clusters 
of galaxies, so smaller errors are expected from a galaxy 
sample. 

The basic method used in this paper is to model the 
observed WMAP maps in all frequencies simultaneously 
using a combination of foreground templates to represent 
the kSZ effect as well as the thermal Sunyaev-Zeldovich 
(tSZ) effect. We also model the contamination from radio 
point sources that may be present in the galaxies, as well 
as Galactic foregrounds. We then analytically marginalise 
all results with respect to primary CMB fluctuations, and 
additionally marginalise over independent monopoles and 
dipoles in each channel. 

In Section [2j we describe the statistical method used 
to fit the templates. In Section [S] we describe the data 
and models used to build the templates for the thermal 
and kinetic Sunyaev-Zel'dovich effects and point source 
contamination. In Section |4j we present and discuss our 
measurement of the kSZ effect. We compare the kSZ bulk 
flow with results using other methods and measure the 
baryon fraction in free electrons in Section[5] In Section[6] 
we discuss future improvements and forecast the errors on 
the components of the kSZ bulk flow from the upcoming 



Planck temperature maps. Finally, Section [7] concludes 
the paper. 



2 STATISTICAL METHOD 

Here, we adopt a full-sky, multi-channel statistical model 
for WMAP7 data. We can model the observed sky tem- 
perature as a combination of multiple physical effects us- 
ing a template based method. The templates that we con- 
sider in this work will model the two Sunyaev-Zel'dovich 
effects, the Point Source contaminations, the Galactic 
foreground emission, the residual temperature monopoles 
and dipoles. Generically, the expected observed signal dc, 
in the channel c, may be written as: 



dc ^ Be i s + Q»,CT.(c)ii,c -I- n 



(1) 



where Be is the beam of the instrument corresponding to 
channel c, n is an instrumental noise component, s is the 
Cosmic Microwave Background (CMB) signal, ai,p are 
scalars and ti^c is the template of signal i in the channel 
c (without beam). (Ti(c) is a projection operator for the 
physical signal i. In case the frequency dependence is not 
known, the projector is the identity. In the other case, the 
projector is a constant. There may be other intermediate 
cases if some channels are at same frequency. 

We assume that the primary CMB signal is a real- 
isation of a Gaussian random field, and that the instru- 
mental noise is also Gaussian. After marginalisation over 
the primary CMB fluctuations, we may write the total 
log-likelihood as (see Appendix [A| : 

X^({Qi,p}) = 

^(dc — ^ a,;,o..(c)-Bcii,c)^C'^_c/(dc'— ^ Oii,CTi(c')Bc'ti^c'), 



c,c' I 



(2) 



with 
and 



D = l + S 



1/2 



"NZ'B,. s 



(4) 



In the above, c and c' run over the available channels of 
WMAP, Nc is the instrumental noise covariance matrix 
in the channel c. The inverse covariance matrix does 
not use the noise covariance matrix (only its inverse) and 
the angular spectrum (and not its inverse). Thus, we do 
not have to worry about how to regulate the inverse of 
these two operators. In the above, we have used pseudo- 
inversion and we have set 

N-^ = MN-'M, (5) 

with M the pixel masking operator. This form has al- 



ready been used in |Wandelt et alT (20041 in the context 
of Gibbs sampling of the CMB fluctuations. 

If we assume that the angular power spectrum of the 
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primary CMB is kept fixed, tlie equation ([2| leads to an 
analytic solution of the a^^pi 



with 
A 



E 



(7) 



c,c 

0-i(c) = p,(Tj(c') = I^ 



Computing ^(i p)_y ,^) requires applying O (A'^t x TVc) 
times the operator . In this paper, we consider ten 
templates, and for one map per channel, three input 
maps. This results in applying C^^ thirty times. We 
do not require more because, for a fixed primary CMB 
spectrum, we may precompute a weighted template i^.c, 
which is sufficient to compute all the terms of the covari- 
ance matrix A and the maximum likelihood estimate a. 
The procedure is detailed in Appendix [B] 

For the practical purpose of the numerical imple- 
mentation, we have separated the fitting of the Galactic 
foreground templates from the fluctuations due to kSZ, 
tSZ and point sources contamination. Doing so, the ma- 
trix A, with 21 X 21 elements, stays the same for all of 
our experiments on the SZ components and point sources. 
We have separated the fitting using Gibbs sampling on 
a set of parameters a^^p assuming the other set is fixed. 
Gibbs sampling provides samples of the posterior distri- 
butions. On the positive side, by selecting samples, we 
have the opportunity of imposing non-Gaussian prior on 
the {ui^p}. On the negative side, we are required to run 
a sampling chain. If the templates are correlated, this 
chain may take time to converge. In practice though, for 
the cases of Section |4j the convergence is achieved in a 
few steps. 



3 TEMPLATE GENERATION 

In this section, we describe the physical models used to 
generate the templates of the three signals considered in 
this work: kinetic Sunyaev-Zel'dovich (kSZ), the thermal 
Sunyaev-Zel'dovich (tSZ) and radio point source contam- 
ination (PSC). 

The signal from the kinetic Sunyaev-Zel'dovich effect 
is our primary goal. This is derived from the galaxies be- 



longing to 2M-I— I- galaxy compilation ( Lavaux & Hudson 



20111, which in turn is based on the 2MASS eXtended 



Sources Catalog ( Skrutskie et al. 2006 1 and publicly avail- 
able redshifts. 

The thermal Sunyaev-Zel'dovich (tSZ) signal is pro- 
duced mostly by clusters of galaxies. As clusters and 
galaxies are correlated, it is necessary to also model the 
tSZ signal and marginalise over it. We model the signal 
using the position and luminosity of clusters in the RBC 
( Kocevski fc Ebeling|20d6 1 catalogue. 

Finally, it may be expected that a subset of local 
galaxies host radio or submillimetre sources that con- 
taminate WMAP data as unresolved point sources. The 
population of galaxies and of unresolved point sources 
may thus be correlated and must also be marginalised 



over. The contamination by unresolved point sources is 
obtained using the information contained in the 2MASS 
extended Sources Catalog (Skrutskie et al. 2006 here- 
after 2MASS-XSC). 

For fitting the contamination of the temperature 
fiuctuations due to Milky-Way emission, we have used 
the templates provided by the WMAP7 collaboration on 
Lambda]^ We have fitted the templates at the full reso- 
lution of TVsidc = 512. 

We now describe the templates in detail. 



3.1 kSZ model 

We assume there is a massive corona of ionised plasma 



around galaxies. Fukugita & Peebles (20061 proposed an 



isothermal model for this corona, which we further mod- 
ify by introducing a hole at the centre of the galaxy. There 
are at least two reasons to do this. First, we reduce the 
possibility of correlating the kSZ signal with radio sources 
located at the centres of galaxies. Second, we expect that 
the plasma close to the centre of the coronae would cool 
and collapse to form the cold gas and stars in galaxies. 
Thus we adopt the following model; 



1 



Pair) 



Qm 2ttG r-2 







if -Rg < r < 7?v, 
otherwise. 



(8) 



with Rv the maximum extension of our gas profile, and 
Rg = 15h~^ kpc. In doing so, we do not attempt to fit 
the part of the galaxy dominated by stars and may avoid 
some point source contamination. The density of elec- 
trons is determined from the gas density, and the Helium 
abundance Y 



ne(r-) = 



Y 



Pair), 



(9) 



with Trip the mass of the proton, and Pg the total mass 



density of the gas, and Y — 0.245 (e.g. Komatsu et al 



20111 



The kSZ effect depends on the electron density and 
the velocity of the electrons with respect to the primary 
CMB, i.e. it is proportional to the momentum. The ve- 
locity of a galaxy in the 2M-I--I- sample can be modelled 
using linear perturbation theory applied to the density 
field. Here instead, for simplicity, we model the fiow field 
to lowest order, i.e. by a dipole or "bulk fiow" with con- 
stant Cartesian components Vi. The average kSZ temper- 
ature change due to a bulk fiow is given by: 



Arksz(n) = — 2rcMB err 



I 



V(«v)2-r2(ft) 



(10) 



with Tc is given by Equation (12 1. After integration, we 



The templates are located at http://lainbda.gsfc.nasa. 



|gov/product/map/dr4/templates_get . cfm . 
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Figure 1. A map of the expected kSZ signal assuming each galaxy has an inward peculiar velocity of 100 km , before 
smoothing with the beam. Only galaxies with cz < 20,000 km have been included for generating this template. The grid lines 
are separated by 30° in the longitude and the latitude directions. 



get: 

ATkszCn) 



tan 



TCMB OrfeO 1 

TrGfieiripC Da 

3 



m 



(11) 

with Q the Heaviside step function, /e being the frac- 
tion of free electrons, fie ~ 2/(2 — F) ~ 1.14, and 
a = 160 km s"'^ the typical velocity dispersion of galactic 
haloes, rc is the projected distance from the centre of the 
galaxy: 

rc{n) = Da{zc) sin [cos~"'"(n ■ he)] = 

DAiz,)y/l-{h-h,)^ (12) 

where Da(2c) is the angular distance of the galaxy at 
redshift Zc, and he is the unit vector in the direction of 
the galaxy. We set the kSZ template to zero at small 
radii rc(n) < -Rg. Finally, we scale -Rv according to the 
luminosity as: 



with L* 



50Ho 



1/3 



(220/i"^ kpc) m-v 



L, 



1/3 



(13) 



3.8 X 10^°/i~^ L0 as given for the best 



Schechter function fit in Lavaux & Hudson (20111. We 



take the normalisation from Fukugita & Peebles ( 2006 1 



The parameter mv is introduced to increase the effec- 
tive virial radius to take into account small halos that 
are not accessible using the 2M+-I- catalogue. We adopt 
my = 2.4 as a default (see Appendix[C]and |Masaki et al.| 



(2011 1) but consider other possibilities in Section 4.2 



We generate three kSZ templates, one for each bulk 
flow Cartesian axis, and then fit simultaneously all com- 
ponents. In Figure [1] we show the kSZ template obtained 
before multiplying by a dipolar bulk-flow and normalised 
to 100 km s~^, in the inward direction. Equivalcntly, it 
would be the kSZ signal that would be produced if all 
selected galaxies were moving at the same speed towards 
us. 

The kSZ template exhibits smaller structures than 
the tSZ signal. Additionally, the kSZ signal does not de- 
pend on the frequency. These two features make the kSZ 
signal far more difficult to detect in the CMB data than 
the tSZ signal. Nevertheless, we note that the kSZ signal 
should dominate the tSZ signal on the scale of the corona 
of a galaxy. The ratio -Rtk between the tSZ signal and 
the kSZ is given by 



-Rtk = 



ATtsz _ 2fcTe 



rUeVc 



(14) 



with V the global velocity of the galaxy, Te the mean 
temperature of electrons in the corona, rrig the mass of 
electrons. With = W*^ K, m^c^ = 511 keV and V = 
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Figure 2. A map of the template for the tSZ signal in the W channel of WMAP, before smoothing with the beam. We use a 
resolution of A^aide = 1024. 



300 km s^^, we get -Rtk ~ 0.34. Additionally, the bulk 
flow gives a dipolar modulation of the galaxy distribution 
on the sky, whereas the tSZ signal would simply trace 
cluster (or galaxy) distribution. Consequently, we do not 
expect a systematic contamination by tSZ at the galaxy 
level. 



3.2 tSZ model 

The tSZ signal is generated by hot and dense ionised gas. 
These conditions are found in the intracluster medium, 
and clusters of galaxies have already been detected using 
this effect (for a review see e.g. 



the model given in the appendix of Afshordi et al. ( 2005 \ 



Birkinshaw 1999 1. We use 



for the temperature and gas density profiles in clusters. 
Given the X-ray luminosity, the RBC catalogue provides 
the virial mass Mvir from which we derive the virial ra- 
dius i?vir = (2GMvir/(Avir^^2))^/^ with Avir = 200, G 

the gravitational constant and H the Hubble constant. 
Using this mass estimate, the model predicts the two pro- 
files assuming hydrostatic equilibrium. The temperature 
change due to the presence of a cluster in direction is 
given by: 



L 



/STtsziy, n) = — 2rcMB cTT/tsz 



I hv 



ds, 
(15) 



with ne(s) and Te(s), respectively denoting the density 
and temperature of the free intracluster electrons at a 
distance s from the centre of the cluster, f the frequency 
at which the temperature is observed. The dependence 
ftsz{x) on the frequency x = hiz/kTcMB is 



ftsz{x) 



(X 
— 
2 



(16) 



rc is the projected distance from the centre of the clus- 
ter. It has the same expression as for galaxies, given in 
Equation ( |12[ ). Note that we stop the integration at twice 
the virial radius of the cluster. 

An example of the tSZ template in the W channel 
of WMAP is given in Figure [2] We have selected all 781 
clusters with redshift z < 0.3, from RBC. With our nor- 
malisation the resulting signal is typically at a level of a 
fraction of a milli-Kelvin. 



3.3 Point source model 



We follow Afshordi et al. (20041 and Jarosik et al. (2011 1 



in modelling the effective temperature of point sources in 
the CMB sky. We remind the reader of their main result: 



s[ni,Pj) 



(17) 



with Fc is a conversion factor between physical flux and 
antenna temperature that depends on the channel c, pj 
is the direction of the j-th pixel, and fii is the direction of 
observation of the galaxy i. We also use 5s{n,m), which 
is one if n and m corresponds to the same pixel of the 
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CMB map, and zero otherwise. a{i'c) accounts for the 
difference between antenna temperature and thermody- 
namic temperature ( Bennett et al.|2003 l: 

4sinh^(a;/2) 

i'rad,c is the average point source luminosity in the chan- 
nel c. A/ps is a normalisation factor to account both for 
pixelisation and the normalisation of the beam. L is the 
luminosity density of galaxies in the considered catalogue. 
We us e the typical value derived by |Lavaux fc Hudson| 
(20111 for normalising: L = 4 x 10* Lq. cf>i is the 



observed flux of the galaxy i: 



(19) 



with Mx,0 the intrinsic luminosity of the sun in the K 
band of 2MASS: Mk,@ = 3.29. [Page etaT] ([20031) reports 



the value for r]f° in the case of free-free emission (scale 
free spectrum with a slope a = —0.1). The point source 
template includes only galaxies with apparent magnitude 
K2Q < 13.5 — Ak (with Ak the extinction derived from 



the reddening (Schlegel et al. 1998 1). The relation be- 



tween the reddening and the extinction is 



Ak = 0.35E: 



(20) 



where the constant of proportionality is obtained from 
the relation between extinction in K band and in V band 
( Cardelh et al.|1989 l. 



3.4 Effective mask 

Since WMAP data are contaminated by strong radio 
sources and foreground emission, it is necessary to mask 
parts of the sky that may contaminate the fitting proce- 
dure. We have used the "Extended Temperature Analy- 
sis" mask (dubbed KQ75 by the WMAP collaboration) 
to remove Galactic contamination. 

In addition, we need to allow for the fact that, due to 
Galactic extinction, galaxies are not uniformly detected 
on the sky. To ensure the completeness of the 2MASS- 
XSC sample, and thus of our template, we also require 
that the Galactic extinction Ak < 0.1. 

The 2M-I— I- sample has more uniform sampling at 
low redshift for Galactic latitudes |6| > 10° because at 
lower latitudes only the galaxies K < 11.5 are properly 
sampled. We thus add two further conditions: Ak < 0.1 
and |ti| > 15° to produce the effective mask. We use 
a stricter than needed to remove also as much galactic 
foreground as possible from the fitting. The final mask 
is shown in Figure |3] The M operator in Equation ([5| 
corresponds to multiplying a pixelated temperature data 
by this mask. 



4 RESULTS 

In this Section, we present our results of the measure- 
ment of the kSZ signal using WMAP7 data. We present 
the significance and the possible shortcomings of the mea- 
surement in Secti0nr4.ll 




Figure 3. Mollweide projection in Galactic coordinates of 
one of the masks used in the analysis. This particular mask is 
generated from the KQ75 WMAP mask, the constraint that 
Ak < 0.1 and |6| > 15°. 




Figure 4. A Mollweide projection in Galactic coordinates 
of the posterior probability of the kSZ bulk-flow direction. 
The dark (light) areas show the direction of bulk flow at 68% 
(95%) probability. We have only shown the constraints corre- 
sponding to the measurement with tSZ signal masked and an 
isothermal density profile (line 1-3 of TableQ. The blue, green 
and red areas correspond to a 50h~^ Mpc, 100h~^ Mpc and 
200/i~^ Mpc respectively, fitting only higher galactic latitudes 
\b\ > 15°. 



4.1 Results from 2M++ 

Our "default" solution is based on the templates and 
mask described in Section [3] with the amplitude of the 
templates left free. We use the best fit angular power 



spectrum provided by the WMAP collaboration (Ko- 
[matsu et al.|[2011[ ) for a ACDM model obtained from 



WMAP7 data alone (Jarosik et al. 20111. Addition 



ally, we independently fit the residual full sky dipole in 
WMAP7 maps, which could otherwise contaminate our 
fit of the kSZ signal, which is also dipolar, albeit localised 
at the positions of galaxies. We produce a set of kSZ tem- 
plates with the 2M+-I- galaxy sample limited at different 
distances from the observer. All templates discard galax- 
ies within 5h~^ Mpc which is the minimum distance at 
which redshift information can be used as an average in- 
dicator of distance. We produce a subsample limited to 
50h~^ Mpc, lOO/i"^ Mpc, 200/i"^ Mpc and fit each galaxy 
sample independently. 

Our results are summarised in Figure |4] Table [l] and 
Table [2) imposing that the Galactic foreground correc- 
tions are only positive. Our results are consistent the 
presence of kSZ signal at the 90% level at line 1. The 
corresponding kSZ map is given in Figure [6] 



© 0000 RAS, MNRAS 000, 000-000 



Bulk flows using the CMB 7 



Table 1. Kinetic Sunyaev-Zcl'dovich estimates from WMAP7 maps. 



Line 


Sky cut 

m >) 


kSZ depth 
(h-i Mpc) 




kSZ 

(km s^-*^) 

Vy 


V^ 




PO 

(%) 


1 


15° 


50 


559 ± 290 


-412 ± 285 


-93 ± 200 


6.1 


11 


2 




100 


305 ± 223 


7 ±222 


-119 ± 151 


2.5 


47 


3 




200 


365 ± 207 


183 ± 206 


-141 ± 143 


5.0 


17 



Note: For all the above measurements, we have used the basic mask of Section |3.4| the tSZ effect as predicted according to our 
pre-beamed tSZ templates, po is the probability of the null hypothesis, namely that there is no bulk flow. The reported has 
three degrees of freedom. Empty field means that the value is kept the same compared to the previous line. For all measurements 
we find consistent tSZ measurement with an amplitude of 1.7 it 0.3, though the bulk of the tSZ pixels is masked out. Without 
masking the amplitude of the tSZ is 1.4 it 0.1. The point source contamination signals are, in unit of 10~^ Jansky, 6 it 2 (Q 
channel), 12 it 7 (V channel) and 40 it 30 (W channel). 



Table 2. Robustness tests and investigation of model systematics on the kSZ signal detection. 



Line Experiment kSZ PO 



(kms-l) (%) 









Vy 








1 


tSZ masked 


559 it 290 


-412 it 285 


-93 it 200 


6.1 


11 


2 


tSZ not masked 


487 it 291 


-218 it 252 


-194 it 190 


5.4 


14 


3 


NEW, tSZ not masked 


555 it 273 


-160 it 262 


-156 it 180 


5.1 


16 


4 


tSZ not masked, my = 10 


573 it 253 


-274 it 243 


-127 it 174 


7.0 


7 


5 


tSZ not masked, cr scaling, my = 2.4 


582 it 386 


-269 it 383 


178 it 302 


3.8 


29 


6 


tSZ masked, a scaling, limit to M» + 1, my = 3.59 


697 it 460 


-483 it 442 


506 it 345 


7.3 


6 


7 


tSZ extended masking, cr scaling, limit to Af* + 1, my = 3.59 


796 it 478 


-407 it 451 


429 it 352 


6.2 


10 


8 


tSZ masked, |fe| > 10° 


531 it 273 


-405 it 259 


-122 it 194 


6.6 


9 



Note: In all the above, we have used a mask based on the KQ75 sky cut, as indicated in Section |3.4| The notation "tSZ masked" 
indicates that we have masked out pixels which were expected to have un-beamed tSZ signal from the cluster sample. "tSZ 
extended masking" indicates that all pixels with a contamination bigger than O.l^K, after beaming in the Q channel, have been 
masked out. All kSZ templates have been generated from the 2MH — \- limited to a depth of 50h~^ Mpc. 



4.2 Tests of potential systematics 

As a check on the foreground template fitting, we report 
m Table |3] our best fit parameters for the Galactic fore- 
grounds (without forcing the parameter to be positive). 
We also show in Table[3]the value of the parameters found 
by the WMAP collaboration, which we have averaged by 
frequency channels. There is a good agreement between 
the two sets of parameters when we include all pixels at 
\b\ > 10°. However, there is some tension on the parame- 
ters of the Tk — Tka when only pixels with > 30° are 
included, but we note that the WMAP collaboration fit- 
ted their models after degrading their map at A'^sidc = 32, 
using all the pixels of the maps. The cut at jbj > 10° is 
thus more consistent with their methodology. We con- 
clude that our foreground fitting at full resolution is con- 
sistent with the findings of the WMAP collaboration. 

After subtracting the Galactic foreground compo- 
nents, we are left with a signal for point-sources which 
is consistent with zero for the V and W channels and 
marginally above zero for the Q channel. The tSZ tem- 
plate has an amplitude of 1.40±0.11, significant at ~ 13(j. 
If the physical model for the tSZ effect was exactly right, 
we should have obtained a value of one. The disagree- 
ment may also come from systematic uncertainties in the 
relation between the X-ray luminosity and the mass of 
the individual clusters. The relation used by |Kocevski fc| 
EbelingI pOOe]) is given by [Allen et"aL] (|2003[) and has an 



uncertainty of ~ 0.20 dex in the zero point. This trans- 
lates to a systematic uncertainty of 50% on the mass of 
clusters derived from the luminosity, which itself leads to 
a ~ 50% normalisation uncertainty of the tSZ template 
map. So, including this uncertainty, our measurements 
are consistent with a tSZ template amplitude of one. 

In Table [2] we investigate how different cuts or mod- 
els affect the results of line 1 of Table [T] The first line of 
this table is a copy of the first line of Table [l] We check 
the impact of pixel masking in line 2. For this experiment, 
we use all the pixels of the map and we assume that our 
model for the tSZ signal is correct in average. This oper- 
ation introduces a marginal change in the kSZ signal. In 
lines 3 to 6, we investigate the impact of the choice of the 
model for the density profile of the plasma corona. Us- 



ing a Navarro- Frenk- White profile ( Navarro et al.||1996 l 
instead of an isothermal profile barely changes the result 
(line 3 vs line 1). We have attempted to scale up the 
virial radius by changing the multiplicative factor my to 
the extreme value my = 10. Increasing this factor may 
be useful for including the presence of smaller haloes that 
are clustered with the galaxy. This operation does not af- 
fect mean and slightly reduces error bars by introducing 
more pixels in the correlation. 

In lines 5 to 7, we have introduced the scaling of 
the velocity dispersions with the luminosity of the galaxy 
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Figure 5. In the right panel, we show the change in direction of the flow resulting from a jackknife test. We removed the 
indicated root pixel of the Healpix mesh for each test. We have used a low-Galactic latitude cut of the mask with |fe| > 10°. The 
kSZ template has been generated with galaxies within a 5-200/i~^ Mpc from the Local Group. The line iso-density contour gives 
the area corresponding to 68% of probability. The root pixels of the Healpix mesh arc given in the right panel. 



Table 3. Foreground contamination of the map estimated from the data. 



Origin 




Dust map 






Tk - Tka map 






Ha map 
ifiK R-i) 






Q 


V 


W 


Q 


V 


W 


Q 


V 


W 


\b\ > 10° 


0.1 ± 0.2 


0.3 ± 0.2 


1.0 ±0.2 


0.37 ± 0.02 


0.18 ±0.02 


0.11 ± 0.02 


1.5 ±0.4 


1.0 ±0.4 


0.7 ± 0.4 


\b\ > 30° 


0.4 ± 0.3 


0.6 ± 0.3 


1.4 ±0.3 


0.26 ±0.02 


0.06 ± 0.02 


0.00 ±0.02 


2.5 ±0.5 


1.8 ±0.5 


1.4 ± 0.5 


WMAP7 


0.20 


0.47 


1.28 


0.23 


0.047 


0.00 


1.22 


0.78 


0.43 



Note: We show the best fit amplitude for each of the templates provided by the WMAP collaboration | |Gold et al.||20lT| l, at a 
resolution of Afgide = 512. The first two lines has been obtained in this work, using the indicated sky cuts. The last line has been 
derived from table 2 of [Gold et ar] | |2011| |. 




Figure 6. kSZ signal corresponding to the mean bulk-flow corresponding to the 50h ^ Mpc cut in the 2MH — [-, limited to |fe| > 15° 
(line 1 of Table [!}. 



© 0000 RAS, MNRAS 000, 000-000 



Bulk flows using the CMB 9 




Figure 7. Test of contamination of the kSZ signal by the 
Galactic plane. We show the direction of our found kSZ bulk 
flow in two cases. We only keep pixels with |6| > 10° {\b\ > 
30° respectively) for the red (blue respectively) contours. In 
both cases, we have limited our sample for generating the kSZ 
templates to 200/1"^ Mpc. 



according to 



a{L) = 160 km s~ 



1/3 



(21) 



This choice effectively enforces that the mass is pro- 
portional to the luminosity, whereas for the fiducial 
choice the mass is proportional to L^/"^, with L« — 
3.9 10^°/i~'^ L0. The first impact of such a choice is just 
to increase error bars (line 5), which is expected as the 
kSZ signal for small galaxies is strongly reduced. Limit- 
ing the sample to galaxies which have an absolute mag- 
nitude lower than M, -I- 1 and adjusting my accordingly 
leads to a strong change in the flow amplitude along the 
Z-direction (line 6). The cause of this change is not un- 
derstood. In line 7, we have tried to remove as much tSZ 
signal as possible by masking out pixels which were con- 
taminated at a level larger than 0.1 fj,K after beaming 
according to our model and the RBC cluster catalogue. 
The kSZ signal does not change compared to line 6. 

Finally, in line 8, we have reduced the masked re- 
gions to include more temperature measurements at lower 
galactic latitudes. Compared to the fiducial experiment in 
line 1, the error bars are slightly smaller in all directions. 
The mean values per directions are in perfect agreement 
with the values in line 1. 

Finally, in the left panel of Figure [5] we report the 
result of a jackknife test obtained by removing the part of 
the sky that corresponds to the indicated root pixel of the 
Healpix mesh. These root pixels are given in the right 
panel of Figure [5] This test shows that, at a given depth, 
the kSZ signal is stable. We note that for deeper measure- 
ments, the signal seems to move along the X-direction 
and thus could be interpreted as a sign of contamination 
by th e WMAP "haze" ( |Finkbeiner|2004[ [Pietrobon et af 
2011 1. Correlations between galaxies and Milky- Way ex- 



tinction properties are likely in the bulge and that could 
correlate the kSZ signal with the haze. We note that the 
effect of this haze would go in the opposite direction as 
what we obtain for the kSZ in all slices: it creates a 
positive temperature fluctuation which would push the 
kSZ signal towards negative X. In our case, the signal is 
pushed towards positive X. Additionally, the push would 
be affected by removal of the the galactic bulge part in 
the Jackknife test, which is not the case. 



4.3 Results from 2MASS-XSC 

Finally, we have estimated the amplitude of the bulk flow 
using the entire 2MASS-XSC limited to < 14.0 as the 
kSZ template. After removing nearby large galaxies and 
parts of the sky that could suffer strong reddening, we 
have 959 554 galaxies with < 14.0. As the 2MASS- 
XSC does not have redshifts, we randomly generated dis- 
tances for each galaxy based on their apparent magnitude 
and an assumed Schechter function for the distribution of 
absolute magnitudes of galaxies. Speciflcally, we used the 
Schechter parameters given in [Lavaux fc Hudson] ( |2011[ ) 
for the 2M-I-I-: a = -0.73, M. = -23.17. The result 
is not strongly dependent on these parameters. We esti- 
mated the kSZ bulk flow from an ensemble of ten reali- 
sations of those distances. The scatter in bulk flow from 
one realisation to the other is ~ 60 km s~^ by compo- 
nent. The median depth of such a survey is 470/i~^ Mpc, 
with a mean distance of ~550/i~^ Mpc because of the 
long tail of the distribution. The results are given in Ta- 
ble [4] and are consistent with no bulk flow on these very 
large scales, with a 95% upper limit of 520 km s~^, after 
correction for random error biasing. The standard devia- 
tion per component of the kSZ signal is ~100 km s~^. We 
have not found any evidence for a bulk fiow that would be 
inconsistent with ACDM. However, our results disagree 
with the bulk flow found by [Kashlinsky et ah] ( |2010[ ) at 
the 4.2(7 level. 



5 BULK FLOWS AND THE FRACTION OF 
FREE ELECTRONS 

While the bulk flow on very large scales is controversial, 
the existence of a bulk flow at depths of < 50 Mpc, 
due to local superclusters such as Virgo, the Great At- 
tractor region and Perseus-Pisces, is not controversial: 
there is good agreement between peculiar velocity sur- 
veys as well as with the predictions of the density held. 
The kSZ signal is the product of the electron density 
times the velocity, i.e. the momentum. Therefore, we can 
take the existence of the bulk flow as a given, and instead 
determine the density of free electrons. 

In order to do so, we must flrst consider the volume 
that is actually being sampled by the WMAP data. Al- 
though our template extends to 200 Mpc, nearby 
galaxy haloes have more weight, and so the bulk flow 
measured is not simply volume limited. We discuss the 
effective volume of the kSZ measurement below before 
discussing the fraction of free electrons. 



5.1 Kinetic Sunyaev-Zel'dovich effect window 
function 

The amplitude and direction of the kSZ effect may be 
interpreted as a bulk flow obtained from the weighted av- 
erage of the peculiar velocities of galaxies. The weighting 
depends both on the selection cut and the limitations of 
the instrument used to observe the CMB. The weighted 
average procedure may be modelled by the convolution of 
the peculiar velocity held by a kernel, called window func- 
tion. We may argue that the bulk flow estimate must have 
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Table 4. Bulk flows at depths > 200 h ^ Mpc from the kinetic Sunyaev-Zel'dovich effect alone. 



Experiment 



Typical depth 
(/i-l Mpc) 



kSZ 
(km s-i) 



v2 




v2 




Xo 


po 


AKAEEK 


Pkaeek 




(%) 




(%) 


4.3 


23 






7.1 


7 


9.1 


2.8 


14.7 


0.2 


16.3 


910-2 


15.2 


0.16 


16.1 


1 10-1 



2MASS-XSC A" < 14 



Kashlinsky et al. ( 2010 i 



470 
175-260 
260-380 
270-530 



160 ± 141 
174 ± 407 
428 ± 375 
352 ± 304 



211 ± 124 
-849 ± 351 
-1029 ± 323 
-713 ± 262 



-11 ± 112 
348 ± 342 
575 ± 316 
652 ± 256 



Note: Xq, and its associated probability po, represent the probability that the flow is zero. XkaeeK' ^'^'^ associated probability 
PkaeeKi is the probability that the results of [Kashlinsky et al. | | [20Tol l are in agreement with our result based on 2MASS-XSC. 




Wq of size R defined as 
Weir) = 



(27r)3/2i?3 



(22) 



Figure 8. We show the window function VV^'^(fc) (compo- 
nent along the X axis of the bulk flow) for our two sample se- 
lections: 50h~^ Mpc (red lines) and 200h~^ Mpc (black lines). 
They have been derived using Eq. ( |D22| of Appendix [d] We 
show in dashed lines the fit to a Gaussian window function to 
the corresponding window functions. The part of the window 
function highlighted by the greyed area may not be trusted 
because of the approximations that we have used. 



a similar window function as tlie one used for maximum 
likelihood analysis of galaxy distance surveys. We may 
assume, e.g., that galaxies have all the same physical size 
~ Rv Their sizes on the sky, in terms of number of pix- 
els, are thus ~ R^/z^ oc with z their redshifts and d 
their angular distances. The number of pixels affected by 
each galaxy in the kSZ template is thus scaling as l/d^. If 
for simplicity we assume that the pixel noise is homoge- 
neous, the weight for the peculiar velocity of each galaxy 
becomes proportional to the inverse squared distance of 
this galaxy. This is exactly the same case for bulk flow 
derived from maximum likelihood estimate from galaxy 
distances. In that case, the peculiar velocity errors scale 
as the distance d, which leads to the same weight per 
galaxy as the kSZ effect. 

We give a quantitative derivation of this window 
function in Appendix |D] where we also discuss the limita- 
tions of this calculation. Figure [S] shows the square of the 
window function of the bulk flow in the Galactic bulge 
direction {X axis). We also show the window function 
corresponding to the convolution with a Gaussian kernel 



The size R is determined such that the second deriva- 
tives of the actual window function and of the Gaussian 
kernel are equal. The measurement, based on galaxies 
within 50h^^ Mpc, seems to correspond to a 18.9h~^ Mpc 
Gaussian window on large scales. Similarly, the measure- 
ment using galaxies within 200h~^ Mpc corresponds to a 
Gaussian window of 3lAh~^ Mpc. We have evaluated the 
window function for the Galactic Y and Z components of 
the bulk flow and we have obtained similar results for the 
probed scales. Because the masked galaxy distribution is 
anisotropic, the window functions are not strictly equal 
for the three components of the bulk flow. The important 
point is that the kSZ signal is dominated by galaxies on 
scales significantly nearer than the limiting depth. 



5.2 Comparison with published bulk flows and 
the free electron fraction 

When only galaxies within 50h~^ Mpc are used in the 
kSZ template, we find \V\ = 533 ±263 km s-^ (after cor- 
rection for "error biasing'Q in the direction I = 324 ±27, 
b — —7 ± 17. This is similar in direction and amplitude 
to recent bulk flows based on classical distance measure- 
ments. For example, at a similar depth, the "Seven Samu- 
rai" (7S) found V = 599 ± 104 km in the dir ection 
Z = 312 ± 11, & = 6 ± 10 ( [Dressier et~aL][T987t using 
the Dn - a distance indicator for early-type galaxies. Ta- 
ble [5] shows three recent bulk flow measurements derived 
using three different methods, along with the character- 
istic scale for averaging the peculiar velocities. We have 
selected these three bulk fiows because they probe simi- 
lar volumes to that probed by the kSZ measurement of 
Line 1 of Ta blefT] For example, th e "SHALLOW" sample 
analysed by Watkins et al. ( 2009[ ) includes measurements 
from surface brightness fiuctuations (Tonry ct al. 2001) 



and the ENEAR sample ( Bernardi et al.|)2b02[ ), an up- 
dated version of the 7S sample quoted above, whereas the 
SFI++ sample contains the latest TuUy-Fisher distances. 



^ The velocity amplitude is a positive number so random er- 
rors will bias the true amplitude upward. We correct for this by 
subtracting off in quadrature the uncertainties in each Carte- 
sian component of the flow. 
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Figure 9. We have drawn the 68% and 95% probability limit 
for the kSZ signal direction and some of the bulk flows reported 
in Tablcjs] Green: Watkins ct al. (2009), limited to 20/i~^ Mpc 
Gaussian. Blue: |Lavaux et al. ( ^2010) , limited to 50/i~^ Mpc 
Top-hat. Red: kSZ signal, template generated with galaxies 
cz < 5000 km s-l. 
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The result from Lavaux et al. (20101 is instead is a pre- 



dicted bulk flow, assuming that 2MASS-selected galaxies 
trace the mass. Figure [9] shows the direction of the bulk 
flow for kSZ measurement (limited to 50 Mpc) and 
for the data sets discussed above. The directions are in 
good agreement. 

Given the agreement in direction, we may assume 
that the amplitude itself is significant. The amplitude 
of the kSZ signal, expressed as bulk flow components 
Vbuik,ksz, may be modelled as the following: 

Vbuik.ksz = /Vbuik, (23) 

with / = /e(cr/160 km s-^)^ /e being the fraction of free 
electrons and <j modelling the unknown normalisation of 
the density profile of the galactic corona. / is taken as 
a free parameter for this subsection, Vbuik the bulk flow 
components as determined by other methods listed in 
Table [s] We may derive the posterior probability of / 
using Bayes identity, which, up to a normalisation, is: 

-P(/|Vbulk,kSZ, Vbuik) oc 
j d'^Vbulk,truc-P(Vbulk,kSzl/, Vbuik, truc)^'(Vbulk, true), 

(24) 

with 

-P(Vbulk,kSzl,f, Vbuik) = 5(,fVbulk,truc, CTkSz), (25) 

and 

P(V 

bulk, true ) = 5(Vbulk, CTbulk). (26) 

In the above expressions, we have used the notation 
0(m, s) for a three dimensional Gaussian probability, 
with the mean rrii and the standard deviation Si for the 
same component, aksz (o"buik respectively) is given by the 
errors on the components quoted in Table [l] (Table |5] re- 
spectively) . Effectively the posterior probability ( |24[ ) may 
be written as a likelihood function 



,f Vi.bull 



V'i.bulk.kSZ 

,2 

i,bulk 

+ log ( /^O-i'.bulk + 



/^<bulk + ^ 

' 72 2 



2 

i,kSZ 



'^'.ksz)] . (27) 



We show t he p osterior distribution of the parame- 
ter / in Figure 10 for different choices of the published 



Figure 10. The posterior distribution for the / factor assum- 
ing different measured bulk flows for a single kSZ measurement 
given by line 1 . The left (right respectively) black solid vertical 
line correspond to / = 1 (/ = Qm/^b = 5.92 respectively). 



bulk flows against one single of our kSZ measurement. 
We are taking the line 1 of Table [l] as the fiducial kSZ 
measurement. Assuming the bulk flow for the SHALLOW 
( |2009[ ), the posterior peaks at 

1.4^ 



of the SFI sample, th i s valu e becomes / 
suming Lavaux et al. (20101 cut at 50^~ 



sample in [Watkins et al.| 

/ = l-itiss (68% confidencelimit). For the bulk fiow 

Mpc top-hat, 

the posterior peaks at / = 1.29lQ g2. ^'^^ information, we 
have shown two extreme choices of choosing the amount 
of baryons in galaxies. The choice / = 1 corresponds to 
considering that galaxies have a dark matter halo and 
all baryons are ionised in the corona. For indicative pur- 



poses, we have shown / = fim.fid/fib = 5.92 in Figure 10 
^^m.fld ~ 0.258. We would measure such a value if the den- 
sity of baryons were equal to the mean density of dark 
matter in the corona of galaxy. We note that all poste- 
riors are in agreement with a choice of / = 1. The ratio 
of the posterior at its maximum to its value for / = in 
the best case (green curve, SHALLOW sample) is 3.97. 
For the blue and red curve, this ratio is 3.93 and 3.08 



respectively. According to Jeffreys ( 1961 1, a ratio of 3.97 



corresponds to substantial evidence against the null hy- 
pothesis. 



6 DISCUSSION AND FUTURE WORK 

There are several ways in which this measurement can be 
improved in the near future. First, improved CMB data 
is anticipated from SPT |S chaffer et al.| ( |2011[ ) and from 
the Planck mission. We discuss below the improvement 
in signal-to-noise expected from Planck. Moreover, we 
anticipate a more detailed treatment of the peculiar ve- 
locity field, an improved understanding of the halo model 
and of the distribution of free electrons around galaxies. 

6.1 Forecast for Planck 

In this section, we estimate the noise covariance matrix 
of the three components of the bulk flows for the maps 
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Table 5. Summary of bulk flows from in the literature with similar depths as the kSZ measurement. 



Reference Scale Fiher Vx Vy \V\ I b 

(h-l Mpc) (km s-i) (km s-i) (km s-i) (km s~i) (°) (°) 



Watkins et al. 


(2009 




Lavaux et al. 


20101 



20 Gaussian 85 ± 27 -228 ±43 53 ± 26 253 ± 40 291 ± 7 12 ± 6 

SHALLOW MLE 190 ± 33 -259 ±29 110 ±23 342 ± 30 306 ± 6 19 ± 4 

SFI++ MLE 117 ± 28 -331 ± 28 96 ± 20 365 ± 27 290 ± 5 15 ± 3 

50 Top-hat 227 ± 77 -270 ± 78 136 ± 78 393 ± 75 310 ± 13 20 ±11 



Note: In the case of "SHALLOW" and "SFI++" , we have used the Maximum Likelihood Estimator (MLE) values quoted in 
IWatkins et"aL] | |2009| . The value given reported as 20h ^ Mpc Gaussian is obtained using the minimum variance estimator applied 
on the "COMPOSITE" sample. 



obtained from the Planck probe. The inverse covariance 
matrix is given by Equation ([7|. Neglecting the impact 
of the other foreground components, assuming the noise 
is homogeneous on the sky and no masking, the covari- 
ance matrix Cb of the kSZ bulk flow becomes diagonal 
in spherical harmonic space and simpHfies to: 



B,,eB,,jCiN-lN;,]^ 



Be', (28) 



with Ti^ij the angular cross power spectrum of the kSZ 
template between the component i and j of the bulk flow. 
For a single channel, this reduces to: 



C 



(29) 



We report our results of the forecast estimate for the 
channel W of WMAP7 in Table [6] In this same table, 
we report the actual errors that we obtained using the 
full method. To make a forecast on the sensitivity that 
we expect using Planck maps, we consider the map at a 
frequency of 100 GIfz. According to [The Planck Collab-| 
(l2006l), 



oration 



the noise at resolution of Ngide = 1024 is 
~ 7 fiK for the temperature. The beam is modelled us- 
ing a Gaussian with a Full- Width-Half-Maximum equal 
to 10 arcmins. We report in the last column of the Ta- 
ble |6] the standard deviation expected for the velocity 
components of the bulk flow with these parameters. 

There is a good agreement between the estimated 
value of the noise for WMAP data and the actual noise 
amplitude that we measure using the complete analy- 
sis of Section [2] This gives us some confidence that the 
projection for Planck is realistic. We note that, mainly 
due to the reduction of the beam width, the noise should 
go down by a factor of two to three by using Planck 
data compared to WMAP data. This estimate does not 
take into account possible additional foreground contam- 
ination. If the bulk flows stay at the same amplitude as 
we have estimated from WMAP data, that could give 
us a ~ 5.4(T detection of the kSZ effect. If the signal 



Component 


WMAP7 


WMAP7 


PLANCK 




obtained 


forecast 


forecast 




(km s^-*-) 


(km s"-*-) 


(km s"-"-) 


X 


290 


305 


120 


y 


285 


277 


111 


z 


200 


224 


93 



Table 6. Standard deviation estimates on bulk flow compo- 
nents for the kSZ templates generated from the 2M++ cata- 
logue cut at a depth 50h~^ Mpc. 



be nonetheless noted that we have neglected all contri- 
butions from the multi-channel analysis in this predic- 
tion. The high-frequency channels would clearly help to 
remove the contribution of the tSZ signal. Additionally, 
the temperature data from the Planck mission have suf- 
ficient angular resolution to probe deeper flows on scales 
up to 200h~^ Mpc scale with better sensitivity. As a final 
note, we assume that the kSZ signal has the same am- 
plitude as what is expected from the bulk flow given by 
the SHALLOW sample in Table [5] In this case, we fore- 
cast that the constraints on the baryon density in galaxy 
coronae from the Planck mission should be A/ ~ 0.2, 
instead of ~ 0.9 as obtained in Section [5.21 



6.2 Model improvements 

There are several ways in which the model for the free 
electron momentum field will be improved. 

(i) The model for the flow field will be updated. In this 
paper, for simplicity, we have modelled the peculiar ve- 
locity field as a simple bulk fiow, a low order moment of 
the peculiar velocity field. However, as suggested by |Ho| 
et al. (20091 and Shao et al. (20111, by using linear per- 



turbation theory and an assumption about the relation- 
ship between mass and light in haloes (linear biasing or 
the halo model), we can predict the peculiar velocities of 
individual galaxies ( Yahil et al.|1986 Strauss et al. 1992| 
Hudson! 1994D to an accuracy of ~ 150 km s~^ ( |Wilhck fc| 
StraussI 19981) or better using a non-linear reconstruction 



( Lavaux et al.|[2b08 20101. These predictions scale with 
the parameter /3 = /(n^)/fe where / ~ fi^'''^ for ACDM 
and fe is a bias factor, and so the kSZ effect will allow us 
to get an independent measurement of this parameter, 
(ii) We are using galaxies as a tracer of the distribu- 
sample), the significance could decrease to 2.2cr. It must tion of free electrons, but we can improve this model in 



goes down to the level of Watkins et al. (20091 (SFI++ 
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several ways. First, the scaling of halo mass and radius 
with galaxy luminosity can be updated using results from 



abundance matching (Marinoni & Hudson 2002 Yang 



et al. 2003 Guo et al 



2010 1 or weak lensing studies of 



low redshift galaxies ( [Mandelbaum et al.||20()6| . Second, 
high-resolution hydrodynamical cosmological simulations 
can be used to get a better handle on the link between 
galaxy halo profiles and their cross-correlations with the 



distribution of free electrons (as in e.g. |Shao et al.||20H 
Peirani et a"L]|2012[ ). 



7 CONCLUSION 

We have developed and applied a method for estimating 
bulk flows from kSZ signals in temperature maps of the 
Cosmic Microwave Background. This method is based on 
a likelihood analysis of the amplitude of a set of template 
maps derived from catalogues of galaxies. Assuming an 
isothermal model for the galaxy gas corona, we have de- 
rived template maps for the kinetic Sunyaev-Zel'dovich 
from the 2M4— I- galaxy catalogue. The thermal Sunyaev- 
Zel'dovich signal is derived in a similar fashion from the 
RBC galaxy cluster catalogue, assuming the model de- 
scribed in |Afshordi et al?] ( |2005[ ) for the density and tem- 
perature of the gas. To estimate the extent of the con- 
tamination by point sources residing in galaxies, we have 
also generated a set of templates from the 2MASS-XSC. 
We have fitted, independently in each channel, the stan- 
dard WMAP7 foreground templates, the monopole and 
the dipole. 

The amplitude of the tSZ signal is measured with a 
high accuracy at 13cr confidence. We have found a 90% 
evidence for the presence of a kSZ signal in WMAP7 data. 
We have not detected any significant contamination by 
point sources. All the results are presented after having 
marginalised over the Galactic foreground emissions. 

Interpreting the kSZ signal as the signature of a 
bulk flow, we derive an amplitude, when only galax- 
ies within 50ft~^ Mpc are included in the template, of 
\V\ = 533 ± 263 km s"\ in the direction Z = 324 ± 27, 
h = —7 ± 17, in excellent agreement with the "Seven 
Samurai" bulk flow (Dressier et al. 1987) and other re- 



cent determinations on similar scales. If we extend selec- 
tion up to 200/i~^ Mpc, we flnd \ V\ = 284 ± 187 km s"^ 
in the direction Z = 26 ± 35, 6 = -17 ± 19. The kSZ 
bulk flow at 50h^^ Mpc is statistically in agreement with 
other measurements obtained from either galaxy peculiar 
velocities, observed or modelled. This comparison can be 
translated into a measurement of the fraction of baryon 
in galaxies. We found substantial evidence of baryons in 
the galaxy corona, compatible with the mean density of 



baryons derived from WMAP7 (Komatsu et al. 20111 



Planck data should be able to reduce the statistical er- 
ror bars on the kSZ measurement by a factor two to three. 
If the signal stays at the present value, the evidence may 
be increased to 5Aa, excluding improvements due to the 
high frequency channels for removing tSZ contamination. 
Pessimistically, if the signal goes down to the value given 



by Watkins et al. (20091, then the evidence would be at 



the 2.2cr level. 
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APPENDIX A: DERIVATION OF THE LIKELIHOOD FROM THE MODEL 

In this Appendix, we derive the Equation ([2| directly from the model assumptions. We start by expressing that the 
primary CMB fluctuations and the instrumental noise by channel are Gaussian. They have known covariance matrices S, 
for the primary fluctuations, and Nc, for the channel c of the instrument. Nc may include masking effect. Nc, as defined 
in Equation ([5| or may both be used in the following. The likelihood of the templates, assuming the amplitudes ai^p 
and the modes s of the primary fluctuations and given the input data, is thus: 



with 



(A2) 



the residual. The above terms can be reordered such that s appears only in one quadratic term: 

X^(s, {a^,p}\{d,}) = {s^ so)^S~^^^DS~'^\s - so) - sls'^^ D'' S'^^so + J2 dl^-^dc (A3) 

c 

with 

and D as defined in Equation (jljl. The likelihood is now split into two parts: 

X^{s, {cti,p}\{dc}) = (s - soY S~'^^^DS~^^'^{s - So) + Y ~ X^Q=i,<^.(c)*i,cj C'cTc' (^c - X^«j,<Tj(c')*j,c' j • (A5) 

The marginalisation according to primary CMB fluctuations consists in dropping the first part, which is the only one 
depending on s. Equation ([2| immediately follows. 

APPENDIX B: ACCELERATING THE PRECOMPUTATION OF THE LIKELIHOOD 

For each template tk,c of the signal k in the frequency channel c, we define a weighted template 

ik,c = s^'^D-''s^'^B,,N;:,hk,c. (bi) 

The Equations ([6| and ([7| may be rewritten in terms of uniquely these weighted templates. We expand the Equation ([3| 
in the Equations Q and Q. After some algebra, we find the expression, numerically simpler, of the covariance matrix 

-^(i,p),(j,"') = 5Z (^icBcNc^B'jj^^, -t^^Bcij^^,^ , (B2) 

c,c' 

0"i(c) = P,CTj(c') = I^ 

and for the maximum likelihood estimate 

Only matrices that are diagonal either in pixel space or in harmonic space are involved in the two above equations. We 
reduce the computational time complexity from 0{Nf x A^^) to 0{Nt x Nc) in terms of the inversion C~^, . 

APPENDIX C: NORMALISATION OF THE GALAXY DENSITY PROFILES 

The mass of a galaxy, for which the profile is limited to the radius R = mvRv, Rv being the virial radius defined in 



Fukugita fc Peebles (20061, is: 



2a 

M(mv) = mv (CI) 

V50GHo 

Assuming an homogeneous distribution of galaxies, with a number of density n*, and that the total amount of dark 
matter is traced by galaxies, the density of dark matter is 

Pm,g = M{m^)nt = — rnv- (C2) 
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Comparing this quantity to the measured mean matter density, 



3Hi 



(C3) 



= (0.41 ± 0.04) mv, 



(C4) 



withn« = (1.13±0.02) lO'^ fe'^ Mpc~^ ( [Lavaux fc Hudson|201l| , fl^ = 0.266±0.028 (WMAP7 alone) and a = 160 km s"^ 
(Fukugita & Peebles 20061. The error bar were added in quadrature to obtain the error on the fraction. This result 



motivates the choice of the effective radius for the matter halo around galaxies. Taking mv = 2.4 yields to an agreement 
between the matter density estimated using galaxies and the one obtained from cosmological probes. 

We have tried to use a Navarro-Frenk- White white profile instead of the isothermal profile. To normalise it in the 
same way, we have enforced that the mass within a virial radius of the isothermal profile is the same as for the isothermal 
profile. 



APPENDIX D: STATISTICS OF THE KSZ ESTIMATOR OF THE BULK FLOW 

For the purpose of understanding the scales probed by the estimator of Equation ||6| in the case of the kSZ effect, we 
model the full observed temperature fiuctuations due to the kSZ signal like: 

dobs(n) = ylksz / dr nl{rn)v{rn) ■ n, (Dl) 
Jo 

with v°'{r) the a component of the velocity field taken at position r, ^ksz the normalisation of the kinetic Sunyaev- 
Zel'dovich, nl{r) the density of electrons at position r and h the direction we are looking in the sky. In this work, we 
extract the bulk flow component by fltting a modulated dipole on the projected electron distribution. Thus the model 
is: 

d„,(V, h) = Aksz dr n:^(rn)V ■ n, (D2) 
Jo 

with Rs the limiting size of the survey used to fit the velocity dipole, and n™(r) is the modelled density of electrons. 

We may write the fitting procedure in spherical harmonic space. In this space, it is possible to marginalise analyt- 
ically according to CMB fluctuation and to take into account beaming effects. The corresponding to the likelihood 
to maximise may be written as: 

1 I I 

X^(V) = ^ — 1%];' - af%(V)| , (D3) 

with Ci — C^^^ + Bj'^Ni, Ni being the noise power spectrum of the instrument, Bi the beam of the instrument, 
the spherical harmonic coefficient of the temperature data dohs(n) 

a-i^ = j A^nYl^{n)dahsin), (D4) 



and the similar relation between af^{y) and cim("-, V). By maximising the likelihood of the bulk flow component, thus 
minimising the above according to the components V", we obtain the following algebraic expression: 

3 

V = Y^{M-^T'''V'' (D5) 

j^a.b = 2^ dV j'^'' dr Yl^{r)n,{rf)f''^ (^J d^s j^' ds y,,™(s)n<,(ss)s*'^ , (D6) 

2^"— / d^fd^s / / drds Yi^{f)Ye.m{s)ne{rr)ne{ss)f''s''v''{ss). (D7) 

^ J Jr=0 Js=0 



with 



V 



We look for the average properties of the estimator given in Eq. ( D5 1 of the bulk flow V. To do this, we want to compute 
the average response of the estimator for different velocity fields and for a fixed positions of tracers. Thus the densities 
rie are kept constant, but the velocity field v(r) is left free. A useful tool to quantify the averaging properties is to 
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consider the covariance {VV*'). Through it is directly related to (VV''): 



d^fd'^sd^r'd^s' 



r=0 J s=Q J J s'=0 



drds yf%(r)Y«,„(s)Yf*,„,(f')Yf/,„/(s') 



X ne(r)ne(s)ne(r')ne(s')f"r' s'^s' {v'^(s)v'^{s')). (D8) 



Assuming that the linear perturbation theory is a good description of the velocity field, that the statistics of density 
fluctuations is purely Gaussian and described by a power spectrum P{k), then the covariance of the velocity field is 



2rr2r2 p+ao 

{V {S)V (S )) = 



dk P{k) / d^fc k^k^e 



2 I. i.'^i.'i ikk-(a-a') 



(D9) 



now rewrite Eq. ( D8 1 



with / — dlogD/dloga, D the growth factor, a the scale factor of the Universe and H the Hubble constant. We may 

(VV") = a^H\f f ^P(fc)VV"'''(fc), 



with 



w^'\k) = iy y f dVd^sdV'd^s' 



67r2 



oo pRs + 



(DIO) 



=0 Jr'^OJa'^O 



drds Yt'^{f)Ye,,n{s)Y/, {r')Yi, ^rn' (s) 



nJr)nJs)nJr')nJs')f''f'''s''s''' / ^^rfcV'-'"-"'^ (Dll) 

./ 4:71 



Up to a multiplication by the inverse of A4'^' , we call this function the window function of the estimator. According 
to Eq. (DIO I, this function is a transfer function to apply on the power spectrum of density fluctuations to obtain the 
covariance of the estimated bulk flow. VV"''' indicates what are the scales probed by the estimator. Other similar window 



functions have been derived in the literature in the context of survey of peculiar velocities of galaxies (Watkins et al. 



20091. Some factorisation may be done: 



W^'ik) = 



4 V V 7T7^fff,mfff',m' / / ne(s)ne(s)s s' / — — k k 



£,c£.£igik-(s-s') 



(D12) 



with 



9e,n 



dr n"Y^%(n)ne(r) 



(D13) 



is exactly the spherical harmonic representation of the kSZ templates that we are using to do the fitting of the 
component a of the bulk fiow in Section [3. 1| The part after the integration symbol in the Equation |D 12| corresponds to 
the complete kSZ signal generated by the whole observable universe. 

The above expression is exact up to the approximation of the linear perturbation theory. To go further we have 
to approximate the behaviour of ne(r). In the following, we will make the following approximation: the galaxies have a 
much smaller size than the considered scales, their distances are large compared to their size, and their density profile 
is going to be simplified to the maximum. We consider that it can be expanded like: 



^galaxy 



ne(r)~ y W^(r-xO, (D14) 



with Xi the position of the i-th galaxy. This sum runs over all galaxy of the observable Universe. W being a function 
with compact support, we approximate it as a small cone: 



Wir~yi,)~AeH{r-x,)e,{f), (D15) 

with Oi{r) = 1 if cos~^(f ■ Xi) < R/xi and zero otherwise, and H{r) = 1 only if \r\ < R, zero otherwise, and A^, the 
mean number density of electrons in the galaxy. This representation of the distribution of electrons in a galaxy is very 
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schematic but sufficient to derive the window function. We now replace the expression ne(r) in Equation (D12I: 



,m V ,m 
j2 



■,mye',m' y ] 



OC /' + 00 



d^sd's' W{s - xOW^(s' - Xj)F,,™(s)y<,,^-(s')s"s''' (D16) 



3d k j^cj^d^ik-(s — s') 



4tv 



p.m. fi' m ' 7. . T 



with i and j running over the galaxies of the observable Universe, and 

K{k) = sinc(fci?)e''="*. 



^d [ 3d^fc 
47r 



(D17) 



(D18) 



(D19) 



We only consider scales much larger than the size of a galaxy, R, thus kR <C 1 and sinc(fci?) ~ 1 to first order in kR. 
We find 



Yt 



47r 



fc"re,(s)ej(s ) 



t,j t,m £' .m' 



3d^fc 

4-7r 



^c^dgik.(xi-Xj) 



XiXj 



2 \ 2 



with 



G"^''(x) = GodxD^"'" + GidxDi^i" 

„ , , „ sin 2: — a; cos x 
Go(x) = 3 ^ 



Gi(x) 



s sin a; + 3a; cos a; — 3 sin x 



and 



— '-■I 



(D20) 
(D21) 

(D22) 

(D23) 
(D24) 

(D25) 
(D26) 



The summation over all galaxies can effectively be truncated to the depth of the survey Rs- In Eq. ( D22 1 the summation 
over the galaxies corresponds to compute the cross correlation between the weighted templates T" and the projected 
density of the Universe of the sky given by the other terms. This correlation is going to be small whenever the galaxy 
i or j is going outside the volume on which we are fitting. We show the amplitude of this effect in Figure |D1[ We note 
that the function is reduced at fc>0.1/i Mpc~^. The shape is preserved for fc^O.l/i Mpc~^. 

We note that the derived window function cannot be trusted on scales k>l/R. The reason lies in the simplistic 
assumption on the density profiles and the approximation on the expansion of the function sine. 

The last missing piece consists in deriving the covariance matrix Al"'''. As we have used an approximation for only 
one part of the derivation of the unweighted window function W"'*', we ought to do the same for the derivation of M"'''' 
even though it is not required. The covariance matrix becomes: 

fRs 



= 47ri? "^"(^^ 

i/xi<Rs 



Ye.m(xi)xi 



(D27) 
(D28) 
(D29) 



with the summation is limited to galaxies with Xi < Rs as indicated by i/xi < Rs- We note that we lose the obvious 
symmetry of A^"''' by making this approximation. Using the expression (D29l, we note that 

^ 3 

W'^k = 0) = ^ M'^-^M"-''. (D30) 

c=l 
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1.0 




Figure Dl. We show the problem of the convergence of the window function for W^'^{k). The fiducial curve corresponds to 
taking the same sample for the galaxies {x^} in Equation | |D22[ | as for the computation of the base template T", in this case 
50h-i Mpc. The deep curve is obtained by considering a sample of galaxies twice deeper for galaxies Xi, thus lOO/i ^ Mpc, than 
for deriving T"'. 



Consequently the weighted window function W°'''{k), defined through the covariance 

3 



{VV'') = a'H^f 



dk 



c,d=l 

is related to W"'^ through: 

c,d=l 

We conclude by acknowledging that W"-'^{k = 0) = (5'''^ 



(D31) 



(D32) 
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